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This  report  discusses,  in  its  first  part,  the  modelling  of  CW  signals  generated  by  an  Acousto- 
Optic  Spectrum  Analyzer  (AOSA).  It  also  shows  how  this  calculation  can  be  related  to  pulse  modu¬ 
lated  signals.  In  its  second  part,  it  discusses  the  processing  of  those  signals  to  detect  and  to  estimate 
the  carrier  frequency  and  the  power  of  an  input  radar  signal.  It  also  proposes  a  system  block  diagram 
to  implement  the  former  two  functions. 


RESUME 

La  premiere  partie  de  ce  rapport  discute  le  calcul  des  signaux  generes  par  un  analyseur  de  spectre 
acousto-optique  pour  des  ondes  continues.  On  montre  aussi  comment  ce  calcul  peut  s’appliquer  a 
des  impulsions  radar.  La  deuxieme  partie  de  ce  rapport  discute  le  traitement  de  ces  signaux  pour 
detecter  la  presence  d’un  signal  et  pour  en  estimer  la  frequence  et  la  puissance.  On  propose  aussi 
une  faqon  d’effectuer  ces  deux  premieres  fonctions. 
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1.0  INTRODUCTION 
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Acousto-Optic  spectrum  analyzers  (AOSA)  are  viable  contenders  to  be  used  as  receivers  for  radar 
electronic  support  measures  (RESM)  applications  especially  because  of  their  high  intercept  probability  and 
their  ability  to  handle  simultaneous  signals.  However,  they  require  a  post-processor  to  extract  the  desired 
information  before  it  is  passed  on  to  a  computer  for  further  processing.  This  report  discusses,  in  the  first 
part,  the  modelling  of  the  output  signal  that  is  generated  by  an  AOSA,  and  in  the  second  part,  the  processing 
of  that  signal  to  detect  the  presence  of  a  radar  signal  and  to  estimate  the  power  and  the  carrier  frequency 
of  that  signal. 


2.0  MODEL 


A  block  diagram  of  the  AOSA  configuration  of  interest  in  this  report  is  shown  in  Figure  1.  The  collimated 
light  wave  generated  by  the  laser  impinges  on  the  Bragg  cell  at  the  Bragg  angle  6g-  The  diffracted  field 
distribution  in  the  frequency  plane  contains  the  frequency  analysis  of  the  input  signal.  Therefore,  an  array 
of  photodetectors  placed  at  the  frequency  plane  can  be  used  to  transform  the  result  of  this  analysis,  which  is 
in  the  form  of  optical  signals,  into  electrical  signals  which  can  be  further  processed  and  analyzed.  The  input 
signal,  which  is  typically  an  electrical  signal  from  an  antenna  feed,  is  transformed  by  the  Bragg  cell  to  an 
acoustic  wave  which  interacts  with  the  optical  beam  with  the  result  that  part  of  the  light  entering  the  cell 
is  diffracted  at  an  angle  proportional  to  the  frequency  of  the  input  signal. 

The  technological  details  and  performance  of  the  various  components  of  this  spectrum  analyzer  are 
beyond  the  scope  of  this  report  but  are  currently  the  object  of  active  research.  The  efficiency  and  bandwidth 
of  the  Bragg  cell,  for  instance,  are  parameters  which  are  sought  to  be  increased,  for  obvious  reasons.  It 
should  be  noted  that  the  input  signal  is  usually  mixed  with  a  local  oscillator,  before  it  is  introduced  to  the 
Bragg  cell,  to  obtain  a  signal  within  the  passband  of  the  transducer. 

Theoretical  formulation  and  experimental  results  of  the  spectrum  analysis  performed  by  this  type  of 
configuration  can  be  found  in  [1],  A  deeper  investigation  of  the  details  of  the  Bragg  cell  diffraction  patterns 
can  be  found  in  [2],  A  mathematical  model  that  is  frequently  used  to  describe  the  signal  processing  of  this 
AOSA  and  the  one  that  will  be  used  throughout  this  report  is  represented  by  the  following  equation  [3][4]: 

Xjt=  /  /  w(z)u(t  -  z)exp(-i2nfz)dz  df  dt  (1) 

where  Xjt  is  the  voltage  produced  at  the  output  of  the  kth  detector  for  the  j ^  time  frame.  In  this  equation, 
u(t)  is  the  input  signal,  te(z)  is  the  window  function  determined  by  the  Bragg  cell  and  the  shaping  of  the 
laser  beam,  /*  is  the  frequency  associated  with  the  k,h  detector,  //(/)  is  a  spectral  weighting  function 
that  describes  the  spatial  response  of  an  individual  detector  element  and  I  is  the  integration  time  of  the 
detectors.  This  equation  implies  that  the  instantaneous  light  intensity  distribution  shining  on  the  array  of 
photodetectors  is  the  instantaneous  magnitude  squared  of  the  Fourier  transform  of  the  signals  propagating 
across  the  Bragg  cell.  Each  photodetector  in  the  focal  plane  spatially  integrates  the  light  intensity  distribution 
and  converts  it  to  currents  which  are  integrated  and  sampled  during  periodic  intervals. 
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2.1  Continuous  Wave  Input  Signal 


For  a  sinusoidal  input,  u(t)  =  Acos(2wf0t  +  <j>),  it  can  be  shown  (see  Appendix  I)  that 


= —n(fk  -  f0)  (2) 

where  H(f)  is  the  convolution  between  the  functions  G  and  H, 

nf)  =  r  H(f  ~  f')G(f')df'  (3) 

J  —  OO 

where  G(f)  is  the  magnitude  squared  of  the  Fourier  transformed  window  function, 


2 


u>(z)  exp(— i2nfz)dz 


(4) 


In  this  report,  it  will  be  assumed  that  H (/)  is  a  rectangular  function  of  unit  amplitude  and  width  B  Hz, 
symmetrical  about  /  =  0.  As  it  was  mentioned  before,  //(/)  is  in  fact  a  spatial  weighting  function  describing 
the  response  of  an  individual  photodetector,  but  there  is  a  linear  one-to-one  correspondence  between  the 
spatial  and  spectral  responses  due  to  the  fact  that  the  angle  of  diffraction  is  proportional  to  the  frequency 
of  the  input  signal.  Although  we  assume  //(/)  to  be  a  rectangular  function,  the  calculations  done  in  this 
report  could  easily  be  redone  for  any  arbitrary  known  //(/). 


2.1.1  General  case 


i 
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From  [5]  we  have  that  the  window  function  u>(z)  can  be  described  by  the  following  equation: 


w(z)  =  exp 


~a(f0)z  -  4 T: 


(f-O  (r-j) 


(5) 


where  a  (in  nepers/sec)  accounts  for  the  acoustic  amplitude  attenuation  and  is  a  function  of  the  input  signal 
frequency,  T  is  the  ratio  of  the  truncated  aperture  over  the  e~2  intensity  width  of  the  Gaussian  profile,  and 
r  (in  seconds)  is  the  transit  time  of  the  acoustic  wave  across  the  Bragg  cell  aperture. 


If  a  =  0  (that  is,  if  we  assume  that  there  are  no  propagation  losses  in  the  Bragg  cell),  we  have  that 
u>(z)  is  a  truncated  Gaussian  profile.  To  see  the  effect  of  the  acoustic  attenuation  on  the  window  function, 
we  can  transform  equation  (5)  into  the  following  form: 
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where  we  can  see  that  the  acoustic  attenuation  causes  a  shift  in  the  peak  position  of  the  Gaussian  profile  as 
well  as  a  decrease  in  the  peak  amplitude.  However,  the  general  shape  of  the  window  function  is  preserved. 

In  general,  there  is  nc  closed  form  solution  to  the  Fourier  transform  of  equation  (5).  So  unless  some 
simplification  or  approximation  is  done,  we  cannot  obtain  a  closed  form  solution  to  G(/).  However,  once  we 
obtain  G(/)  we  can  easily  obtain  'H(f)  by  numerically  integrating  G(f)  over  finite  periods  since  we  assume 
that  //(/)  is  a  rectangular  weighting  function. 

In  Appendix  II  we  show  how  the  rectangular  rule  for  numerical  integration  can  be  used  to  obtain  an 
aliased  version  of  G(f).  The  resulting  normalized  numerical  equation  is: 


Ga(v)  =  r2  {[TZe  (AT(t>))]2  +  [7m  (/r(v))]2}  (7) 

where  v  =  fr  and 

K  El'o1  h m  cos  (V)  -  m  sin  (8) 

where  N  is  the  number  of  points  from  w(z)  that  will  be  used  to  calculate  one  sample  of  G(/).  We  see  that 
the  above  algorithm  is  closely  related  to  the  DFT  algorithm  except  for  the  fact  that  it  can  be  used  to  obtain 
samples  of  G(/)  at  any  frequency  instead  of  obtaining  those  for  a  fixed  set  of  frequencies.  Normalizing  u>(z) 
to  the  following  equation: 


w(x)  =  exp 


-Lx -AT2 


(*-£)  '"‘(—0 


(9) 


where  L  =  or  and  x  —  f,  makes  it  easy  to  calculate  equation  (7).  Ga(v)  is  periodic  with  a  period  of  N 
with  respect  to  v,  so  it  should  only  be  used  to  calculate  samples  of  G(v)  for  the  range  =£-  <  v  <  %.  Since 
Ga(f)  is  an  aliased  version  of  G(/),  N  should  be  made  large  enough  to  make  this  error  insignificant  and  we 
should  not  attempt  to  calculate  samples  of  G(f)  which  are  lower  than  a  certain  value.  The  value  of  N  and 
the  lowest  value  of  G(/)  that  we  attempt  to  calculate  using  Ga(v)  depend  on  G(f),  so  determining  this  may 
require  some  trial  and  error.  When  a  =  0  and  T  =  0,  u/(z)  is  simply  a  rectangular  window  and  G(/)  is  a 
squared  Sine  function.  For  that  case,  the  sidelobe  level  is  down  approximately  50  dB  from  the  peak  at  a 
frequency  of  ^2.  For  all  other  cases,  the  sidelobes  will  decrease  even  more  rapidly. 

Figure  2  shows  the  window  function,  u;(z),  for  the  case  of  or  =  0.5,  T  =  1.  Figure  3  shows  the  function 
G(/)  calculated  from  equation  (7)  using  N  =  200  for  the  same  case  of  a r  =  0.5,  T  =  1.  Finally,  Figure  4 
shows  H(f)  for  that  same  case  assuming  the  photodetectors  have  a  bandwidth  of  B  =  £.  It  should  be  noted 
that  H(f)  is  symmetrical  about  /  =  0  even  though  this  is  not  shown  on  Figure  4.  This  last  Figure  is  the 
result  of  numerically  integrating  the  function  G(/).  For  that  reason  it  is  useful  to  have  a  numerical  algorithm 
that  can  calculate  samples  of  G(/)  at  any  frequency  because  most  numerical  integration  programs  require  a 
function  which  can  do  that. 
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2.1.2  Uutruncated  Gaussian  windowing 


The  preceding  section  showed  how  the  output  of  the  AOSA  could  be  calculated  using  an  algorithm 
closely  related  to  the  DFT  algorithm.  But  it  has  long  been  recognized  that  the  DFT  is  computationally 
expensive.  In  fact  it  is  presumed  that  Carl  Friedrich  Gauss,  the  eminent  German  mathematician,  developed 
an  algorithm  that  could  simplify  its  calculation  as  early  as  the  year  1805  [6]. 

It  is  therefore  worthwhile  to  note  that  if  we  assume  that  a  =  0  and  that  the  Gaussian  shaping  is  not 
truncated  by  the  physical  size  of  the  Bragg  cell,  then  we  can  get  a  closed  form  expression  for  G(/).  Harms 
and  Hummels  have  used  this  result  in  [4]  and  have  calculated  the  transform  of  tu(z)  by  means  of  a  line 
integral  in  the  complex  plane. 


If  we  let 


tu(z)  =  exp  — 4T 


(Hfl 


G(v)=r2-^exp 


/-W\ 
V  2  r2  ) 


where  again  v  =  /r.  This  greatly  simplifies  the  calculation  of  H(f)  but  it  should  only  be  used  when  T  is 
large.  In  [4]  this  approximation  for  T  =  1.63  is  used  and  it  is  claimed  that  the  resultant  error  is  negligible. 

Figure  5  shows  the  window  function,  w(z),  for  the  case  ofT  =  2  while  Figure  6  shows  the  function  G(/) 
calculated  using  equation  (11).  Finally,  Figure  7  shows  7f(/)  for  that  same  case  assuming  the  photodetectors 
have  a  bandwidth  of  B  = 


2.1.3  Truncated  exponential  windowing 


If  we  assume  that  T  =  0,  that  is,  if  we  assume  that  the  laser  has  no  Gaussian  shaping  but  we  still  want 
to  take  into  account  the  fact  that  there  is  attenuation  of  the  acoustic  wave  as  it  propagates  through  the 
Bragg  cell,  then  vu(z)  becomes: 


tu(z)  =  exp(— az)  rect  ^ 


and  it  is  easy  to  show  through  simple  integration  that: 


,[l-2e  L  cos(2?rv)  +  e  2t] 

G(v)  =  t - - - * - - 

L2  +  (2ttv)2 


FIGURE  7  -  H(f)  FOR  ar  =  0,  T  =  2,  tB 


Figure  8  shows  the  window  function,  u/(z),  for  the  case  of  ar  =  0.5  while  Figure  9  shows  the  func¬ 
tion  G(f)  calculated  using  equation  (13).  Finally,  Figure  10  shows  H(f)  for  that  same  case  assuming  the 
photodetectors  have  a  bandwidth  of  B  = 


2.1.4  Rectangular  windowing 


As  it  was  mentioned  earlier,  if  we  ignore  the  acoustic  loss  of  the  Bragg  cell  and  the  Gaussian  shaping 
of  the  laser,  then  w{z)  is  simply  a  rectangular  window  whose  duration  is  determined  by  the  physical  size  of 
the  crystal  or  the  width  of  the  light  wave  impinging  on  the  cell.  For  that  case  it  is  easy  to  show  that  w(z) 
is  as  shown  in  Figure  11  and  G(f)  is  given  by  the  following  equation: 


G(f)  =  r 


2  s«n2(2jr/r/2) 
(2*/r/2)2 


which  is  shown  in  Figure  12.  The  resulting  function  'H.(f)  for  rB  =  1  is  shown  in  Figure  13. 


2.2  Pulse  Modulated  Input  Signal 


The  preceding  section  has  dealt  with  the  case  where  the  input  is  a  sinusoidal  signal  or  what  is  also 
called  a  continuous  wave  (CVV)  input  signal,  in  the  radar  context.  However,  most  radar  signals  are  pulse 
modulated  and  one  might  ask  how  relevant  is  the  preceding  section. 

It  turns  out  that  this  depends  on  the  relationship  between  the  duration  of  the  window  function,  r,  the 
integration  time  of  the  photodetectors,  /,  and  the  pulse  width  of  the  radar  signal.  For  example,  if  the  width 
of  the  radar  pulse  is  several  times  larger  than  r  and  /,  then  the  results  of  the  preceding  section  are  perfectly 
valid  to  describe  the  outputs  from  the  photodetectors  for  most  of  the  integration  time  frames  occurring 
during  the  pulse. 

One  might  ask  about  the  other  extreme  when  the  width  of  the  radar  pulse  is  smaller  than  r.  Lee  [1] 
has  done  some  useful  calculations  in  that  respect.  There  it  is  mentioned  that  there  are  two  general  effects 
of  the  window  function  on  the  time-integrated  illumination  pattern  shining  on  the  photodetectors: 

(1)  a  broadening  effect  on  the  main-lobe  width 

which  inevitably  implies  an  increase  of  the  side-lobe  rolloff,  and 

(2)  a  smoothing  effect  on  the  side-lobe  structure. 

It  should  be  noted  that  the  integration  time  used  in  [1]  was  variable  and  was  equal  to  the  total  duration 
of  time  for  which  the  pulse  is  propagating  through  the  cell  and  interacting  with  the  light  wave.  But  the 
interesting  result  is  that  when  the  width  of  the  pulse  is  equal  to  r/4  or  smaller,  then  the  time-integrated 
illumination  pattern  can  be  very  closely  approximated  by  a  Sine2  function  with  the  first  zero  located  at 
l/PW ,  where  PW  =  pulse  width. 


FIGURE  8  -  WINDOW  w(z)  WHEN  ar  =  0.5,  T 


FIGURE  9  -  G(f)  WHEN  ar  =  0.5,  T 


FIGURE  10  -  H(f)  FOR  or  =  0.5,  T  =  0,  rB 


FIGURE  II  -  WINDOW  u/(r)  WHEN  ar  =  0,T 
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We  can  therefore  get  an  appreciation  of  the  validity  conditions  for  the  results  of  the  previous  section. 
Basically,  the  spectrum  of  the  input  signal  is  unaffected  by  the  window  function,  tu(2),  when  the  width  of 
the  radar  pulse  is  smaller  than  r/4  and  it  is  greatly  affected  by  it,  in  the  way  that  it  was  described  in  the 
previous  section,  when  the  width  of  the  radar  pulse  is  larger  than  r. 


3.0  PROCESSING 


The  preceding  section  has  shown  how  we  can  calculate  H(f)  for  different  system  parameters.  Of  course, 
once  we  have  H(f)  it  becomes  easy  to  calculate  the  outputs  of  the  AOSA  because  these  two  are  related  as 
shown  in  equation  (2).  As  an  example,  the  outputs  of  the  AOSA  are  shown  in  Figures  14  and  15  for  an 
input  signal  u(t)  =  A  cos(2n- fot  +  4>).  Figure  14  shows  the  output  when  tB  -  3/4  for  the  case  where  tu(2)  is 
a  rectangular  function.  Figure  15  shows  the  output  when  tB  =  6  for  that  same  case.  Both  of  these  Figures 
show  the  output  when  f0  corresponds  to  one  ft  as  well  as  when  /0  is  exactly  in  the  middle  of  two  fk’s. 

In  this  section,  we  deal  with  the  processing  of  the  outputs  of  an  AOSA  to  detect  the  presence  of  a  signal 
and  to  estimate  the  power  and  frequency  of  that  input  signal.  We  will  approach  these  three  problems  using 
classical  detection/estimation  theory  as  expounded  by  Van  Trees  [7]. 


3.1  Detection 


If  we  let  R  =  {ri,r2,r3,  ...,rN}  be  the  received  vector  where  the  r/s  are  the  pixel  output  values  for  a 
given  frame,  then  we  can  consider  this  problem  as  a  binary  hypothesis  test: 

Ho  :  r,  =  rii 

Hi  :  r,  =  m,  +  n<,  i  =  1,2,3,  ...,N 

where  the  rij’s  are  assumed  to  be  independent,  zero-mean,  Gaussian  random  variables  with  variance  <y2  and 
the  m/s  are  defined  by 

A2/ 

mi=—Hi  (16) 

where  the  Hi  s  are  the  samples  from  H(f)  which  correspond  to  the  pixel  values  of  our  given  system. 

We  will  assume  for  now  that  the  Hi's  are  known.  That  is,  we  are  seeking  a  test  to  decide  whether  a 
particular  signal  is  present  or  not.  Even  though  we  use  a  notation  which  is  consistent  with  the  previous 
section,  the  following  analysis  is  more  general  and  applies  to  the  detection  of  other  signals  which  have  not 
been  dealt  with  in  that  section. 

In  practice,  the  noise  distributions  will  not  have  a  mean  of  zero  but  will  rather  be  characterized  by  what 
is  often  called  “fixed  pattern  noise”.  However,  as  is  shown  in  [8],  it  is  easy  to  determine  the  means  and 
subtract  them  so  that  this  zero-mean  assumption  is  reasonable. 
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Since  in  the  RESM  context  we  have  no  a  priori  knowledge  of  the  probabilities  for  the  occurrence  of 
hypothesis  H0  or  Hi,  we  cannot  make  use  of  the  Bayes’  criterion.  If  we  let  Pf  be  the  probability  of  a 
false  alarm  (i.e.  the  probability  that  we  say  the  signal  is  present  when  it  is  not)  and  Pd  the  probability  of 
detection  (i.e.  the  probability  that  we  say  the  signal  is  present  when  it  is),  then  an  obvious  criterion  would 
be  to  constrain  Pf  and  design  a  test  that  will  maximize  Pp-  This  is  known  as  the  Neyman-Pearson  criterion 
and  it  is  shown  in  [7]  that  it  is  satisfied  by  the  likelihood  ratio  test 


where  the  likelihood  ratio  is  defined  as 


HR)  =  P\jrr-\  (18) 

p(R  I  Ho) 

where  p(R  \  Hi)  is  the  conditional  probability  of  the  received  vector  R  given  that  hypothesis  Hi  is  true  and 
where  A  is  defined  by  the  following  equation: 


PF  =  J”  p(\(R)  \  H0)dA(R) 


It  is  easy  to  show  that  for  our  problem  as  it  is  defined  in  equation  (15),  the  above  conditional  probabilities 
are  as  follows: 


P(«  I  Wo)  =  n  ^=^5  exp  (^-) 


hence  the  likelihood  function  can  be  written  as 


A  (R)  —  exp 


TlLr?-r.L(ri-mif 


or  even  as 


A  (R)  =  exp 


2<r2 


But  because  the  natural  logarithm  is  a  monotonic  function  and  both  sides  of  equation  (17)  must  be 
positive,  the  likelihood  ratio  test  is  equivalent  to  the  test 


WWW 

which  is  very  useful  given  the  form  of  A(#)  as  shown  in  equation  (23). 

Indeed,  we  can  say  that  a  test  satisfying  the  Neyman-Pearson  criterion  for  our  problem  is  as  follows: 

2E.-=i  r»m«  ~  Ei=i  mi  >  to* 


which  is  equivalent  to  the  test 


W  // 

^rjmi  >  7  (26) 

i=l  H0 

where  7  is  a  constant. 

It  should  be  noted  that  hypothesis  Hi  is  a  composite  hypothesis  because  it  contains  the  “unwanted” 
parameter,  A  and  the  test  of  equation  (26)  cannot  be  used  unless  we  know  the  value  of  A.  However,  since 
m<  =  ,  we  can  reduce  the  above  test  to  the  following: 

5>W«  %  V  (27) 

«=1  Ho 


because  A2 1  is  a  positive  quantity.  If  we  let  Y  =  r,7fj,  then  the  above  test  becomes 


•L *y| 


Ho)dy 


The  above  test  is  what  is  called  a  “uniformly  most  powerful”  (UMP)  test  because  it  is  completely  defined 
(including  threshold)  without  knowledge  of  the  parameter  “A" ,  which  is  the  amplitude  of  the  input  signal. 
Of  course,  the  performance  of  this  test  will  be  a  function  of  this  parameter,  as  will  be  shown  shortly. 

In  order  to  find  the  performance  of  this  test,  we  note  that  if  H 0  is  true,  then  Y  is  a  Gaussian  random 
variable  with  zero  mean  and  variance  <ry  =  <r2  Ti2.  And  similarly,  if  H\  is  true,  then  Y  is  a  Gaussian 
random  variable  with  mean,  Y  =  J2i=i  ™<Hi  =  H2  and  variance,  a\  —  <r2  5Zi=i  From  this 

we  have  that 


PF=r^’xp(.^)d‘ 
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where 


and 


(31) 


(32) 


Pd 


V^aY  CXP 


~-(*-y)a 

2  <Ty 


dx 


(33) 


or 


Pd  =  er/c 


fir) 


(34) 


Or  if  we  use  the  error  function  notation  of  the  IMSL  library  we  have  that: 


Pf  =  j  erfciMSL 


C/D 


(35) 


and 


i  .  /v-y\ 


Pd  =  2 erfcIMSL 


(5 


\Z2(Ty  ) 


(36) 


where 


2  y00 

erfciMsi(t)  =  e\p(—t2)dt 


(37) 


Figure  16  shows  the  performance  of  the  test  of  equation  (27)  as  a  function  of  the  parameter 


x 


<T 


N 


which  serves  as  a  figure  of  merit.  Indeed  x  is  the  S/N  ratio  times  the  integration  time  of  the  photodetectors(  I) 
times  a  factor  which  quantifies  the  performance  gain  due  to  the  matched  filter.  We  can  see  that  the  perfor¬ 
mance  gain  due  to  the  matched  filter  depends  on  the  number  of  pixels  that  are  used  in  the  algorithm. 


a 


One  would  like  to  know  how  the  performance  of  the  test  of  equation  (27)  is  affected  when  the  Hi’s  used 
in  the  test  are  different  from  the  Hi's  that  correspond  to  the  signal.  This  could  happen  when  the  frequency 
of  the  input  signal  is  different  from  the  one  we  assumed  or  when  the  signal  is  actually  a  short  pulse  or  when 
the  window  function,  t/;(z),  that  we  used  to  calculate  the  Hi’s  was  not  accurately  measured  or  when  the 
model  used  to  calculate  the  Hi’s  does  not  accurately  represent  the  processing  of  the  AO  processor  or  for  any 
other  reason. 

To  calculate  this  performance  degradation,  let 


N 


t=l 


Hi 

> 

< 

Ho 


(38) 


be  the  test,  where  the  Hi’s  are  used  in  the  test  but  do  not  actually  correspond  to  the  Hi's  of  the  signal.  In 
this  case,  if  Ho  is  true,  Y  is  a  Gaussian  random  variable  with  zero  mean  and  variance  cr?  = 

But  if  Hi  is  true,  then  Y  is  a  random  variable  with  mean 


N 


A3 1 


N 


E{Y}  =  =  ~^H,H 

»=1  1  =  1 


and  variance  <r?  =  fr2£^i  H3.  It  is  easy  to  show  that  the  performance  of  the  test  of  equation  (38)  is  the 
same  as  the  performance  of  the  UMP  test  as  shown  in  Figure  16  except  that  for  this  case 


x  = 


a’ i pi,  fun. 


and  hence  the  performance  degradation  is  determined  by  the  term  H,H,  which  can  be  seen  as  a  figure 
of  merit  to  evaluate  the  degradation  of  the  proposed  test. 

It  should  be  noted,  in  concluding  this  section,  that  in  the  RESM  context  one  has  no  a  priori  knowledge 
of  the  frequency  of  the  input  signal.  However,  one  could  use  the  proposed  test  to  decide  on  the  presence  of 
many  different  frequencies.  This  could  be  done  by  using  an  FIR  filter  where  the  tap  coefficients  are  matched 
to  the  Hi’s  that  are  expected,  the  implementation  of  which  is  discussed  in  [8]. 


3.2  Power  Estimation 


As  it  was  mentioned  in  section  3.1,  hypothesis  Hi  is  a  composite  hypothesis  because  it  contains  the 
unknown  parameter  “A”,  which  is  the  amplitude  of  the  input  signal  (c.f.  equations  (15)  and  (16)  )  However, 
it  is  possible  to  estimate  the  value  of  the  parameter  “A”  from  the  r,’s  of  the  received  vector  R.  But  since 
the  power  of  the  input  signal  is  usually  of  more  interest  than  the  amplitude,  we  will  estimate  “A2”. 

There  are  many  algorithms  that  could  be  used  to  estimate  “A2”,  however,  one  that  has  an  important 
property  is  the  maximum-likelihood  (ML)  estimate.  Indeed,  it  is  mentioned  in  [9]  that  if  the  maximum- 
likelihood  estimate  is  unbiased,  then  it  will  satisfy  the  Cramer-Rao  bound  with  an  equality.  An  estimator 
which  satisfies  the  Cramer-Rao  bound  with  an  equality  is  known  in  the  literature  as  an  efficient  estimator. 
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Now  to  find  the  ML  estimate  for  A 2  we  must  solve, 


dln{p(R\A 2)} 


Substituting  equation  (21)  in  equation  (39)  and  simplifying  gives  us  that 


Substituting  the  m*’s  gives  us 


-  m ■ )2/2q,z 

8A 7 


ozL&iu-'i)1 


Differentiating  and  simplifying  we  obtain 


£(^’-,k)=o 

•  =  1  v  / 


From  which  we  get  the  desired  ML  estimate: 


-  4  <C»=i  r’^' 

cl  Ml  —  —  _  k. 


1  £?=,*? 


If  we  seek  to  find  the  mean  of  the  ML  estimate,  we  have  that 


e{a\,l } 


which  becomes 


g{^ML}  =  7„,v'-TE^C^» 

1  2^i  =  l  i  =  l 


which  can  easily  be  calculated  to  be 


E{A\,l)  =  a7 


Hence  we  see  that  the  ML  estimate  is  unbiased  and  we  can  therefore  conclude  that  it  is  an  efficient 
estimate  of  the  parameter  A2 .  That  is,  it  is  impossible  to  find  an  other  unbiased  estimate  of  A 7  which  has 
a  smaller  variance  than  the  ML  estimate  of  A7 . 
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Of  course,  we  could  use  our  estimate  of  A 2  to  calculate  the  likelihood  ratio  of  equation  (26)  and  therefore 
obtain  what  is  called  a  generalized  likelihood  ratio  test.  However,  this  would  not  be  a  good  thing  to  do  since 
we  were  able  to  obtain  a  UMP  test  for  this  problem  and  the  generalized  likelihood  ratio  test  is  bound  to 
have  a  poorer  performance. 

However,  for  sake  of  completeness,  it  can  be  shown  that  the  generalized  likelihood  ratio  becomes: 


(Y,)2  f  7j 

Ho 


(47) 


where  y2  is  a  constant  and  Y„  =  J2iLi  ri%-  U  can  further  be  shown  that  the  performance  of  this  test  is 
such  that 


Pd  =  erfc  +  erf  ( 


(48) 


and 


Pf  =  2  erfc 


(£) 


(49) 


where  Y}  is  the  mean  of  Ya  and  a\  is  the  variance  of  Yt.  Figure  16  also  shows  the  performance  of  that  test 
where  again 


x  = 


A2I 


N 


i=i 


as  for  the  UMP  test. 


3.3  Carrier  Frequency  Estimation 


In  section  3.1  it  was  mentioned  that,  in  the  RESM  context,  we  have  no  a  priori  knowledge  of  the 
probabilities  for  the  occurrence  of  hypothesis  Ho  or  //j.  But  even  worst  than  that,  we  also  don’t  have  any 
a  priori  knowledge  of  the  probability  distribution  function  for  the  carrier  frequency  of  the  input  signals. 
However,  if  we  are  given  that  a  signal  is  present,  we  can  express  the  likelihood  function  as 


P(R\fo)-Y[j-ff2exp{{  ,2a3‘)  } 


(50) 


where  the  frequency  dependence  rests  upon  the  rrij's  which  are  assumed  to  be  known  for  different  carrier 
frequencies,  /o-  Or  to  simplify  the  calculations,  we  can  obtain  the  log  likelihood  function  which  is 
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-JVIn(v/&)-^(’’l2gTl)~  (51) 

1=1 

By  definition,  the  ML  estimator  of  the  carrier  frequency  is  the  frequency  that  will  maximize  equation 
(51).  But  since,  in  general,  there  exists  no  closed  form  expression  for  the  m,1  s,  it  is  not  possible  to  get  an 
expression  for  the  ML  estimator.  One  way  to  counter  this  problem  is  to  substitute  the  m, ’ s  which  correspond 
to  all  the  possible  frequencies  in  equation  (51)  and  note  which  frequency  actually  maximized  it.  If  there 
are  an  infinity  of  possible  frequencies,  then  instead  we  can  separate  the  range  of  frequencies  in  a  number  of 
bins  which  form  a  fine  mesh  over  that  range.  In  that  case,  the  number  of  bins  could  be  determined  by  such 
factors  as  the  frequency  resolution  desired  or  the  computing  power  available  or  affordable. 

It  is  clear  that  equation  (51)  is  maximized  when 


\R-  M\2  =  |/?|2  +  |M|2  -  2(R,A/) 


(52) 


is  minimized.  In  general,  | A / 1 2  is  not  constant  for  all  frequencies  but  is  periodic  with  period  B,  which  is 
the  bandwidth  covered  by  one  photodetector.  Hence,  if  we  calculate  equation  (52)  for  a  set  of  frequencies 
which  are  separated  in  frequency  by  multiples  of  B,  then  equation  (52)  will  reach  a  minimum  when  the  term 
(R,H)  reaches  a  maximum.  This  is  an  interesting  result  in  light  of  the  fact  that  this  term  corresponds  to 
the  UMP  detection  test  that  was  obtained  in  section  3.1. 

Indeed,  a  block  diagram  of  a  system  to  detect  the  presence  of  signals  and  estimate  the  carrier  frequency 
of  those  signals  could  be  as  shown  in  Figure  17.  The  matched  filter  could  be  implemented  by  an  FIR  filter 
where  the  tap  coefficients  are  the  calculated  Hi's.  If  many  different  types  of  signals  could  be  present,  we 
could  select  a  filter  which  gives  the  smallest  overall  degradation  by  using  the  result  obtained  in  section  3.1. 
The  output  of  the  matched  filter  can  be  compared  with  a  threshold  comparator  to  decide  if  a  signal  is  present 
or  not.  And  finally,  of  all  the  outputs  of  the  FIR  filter  which  are  above  the  threshold,  the  highest  one  could 
be  used  to  estimate  the  carrier  frequency  of  the  input  signal.  It  should  be  noted  that  the  peak  detector  block 
should  also  have  the  capability  of  recognizing  more  than  one  signal  per  integration  frame  which  could  often 
be  the  case  for  this  application. 


4.0  CONCLUDING  REMARKS 


We  have  discussed  the  modelling  and  processing  of  the  signals  from  an  AOSA.  The  modelling  has  been 
done  for  CW  signals  but  we  have  shown  that  it  can  be  applied  to  pulse  modulated  signals  as  well.  The  case 
of  biphase  coded  signals  has  not  been  addressed.  The  analysis  has  been  kept  general  as  much  as  possible  so 
that  the  results  can  be  tailored  to  specific  systems.  We  have  proposed  the  block  diagram  of  a  system  that 
could  be  used  to  detect  and  estimate  the  carrier  frequency  of  an  input  radar  signal. 
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FIGURE  17  -  DETECTION/ESTIMATION  SYSTEM 
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We  have  defined, 


r(j  +  \)I  roo  rco  2 

Xjt=  /  /  H(f  -  ft)  /  w(P)u(t  -  /?)exp(-i2jr//?)d/?  df  dt  (I-  1) 

J  jl  J  —  OO  J  —  oo 


which  can  be  rewritten 


+  r°°  /-oo  <•» 

A :,*=/  /  /  / 

«/  j  /  J  —  oo  J  —  oo  J  —  oo 


exp  [-«2ff/(o  -/?)]</»  dfidfdt  (1-2) 


Now,  if 


u(<)  =  Acos(2irfot  + 


(/-3) 


u(t  -  a)u*(<  -  /?)  =  A  cos  [2ir/0(<  -  a)  +  <£]  ■  ,4  cos  [2jr/0(t  -  /?)  +  <£] 


(/-4) 


which  can  be  rewritten  as 


^  {cos  [Air fQt  -  2irf0(a  +  0)  +  2<f>\  +  cos  [2 vf0(a  -  /?)]} 


(/-5) 


Hence, 


42  y(i  +  l)/  r°°  r°o  ,oo 

Xjk  -  T  /  /  /  /  ^(/ - /k)^(o)^‘(/?){cos[2T/o(a- ^)] 

^  j  jl  J  —  oo  J  —  oo  J  —OO  ^ 


Jjl  J  —  oo  J  —  oo  J —oo 


+  cos  [4?r/0t  -  2Jr/o(a  +  0)  +  2(/>]  j  exp[-i27r/(a  -  /?)]  dad/ddf  dt  (/  -  6) 


Xjk  =  ^r(  j  f  H(f  -  fk)w(a)w*(0)exp[-i2irf(a  -  /?)] 

^  J  —  oo  J —oo  J  —  oo 

r«+‘V  r  s 

y  |cos[2t/0(o  -/?)]  +  cos  [4x-/o<  -  2irf0(a  +  0)  +  2<j>}jdt  da  d/3df  (I-  7) 


exp  [-t'2x/(a  -  0)]  J  |cos  [2irf0(a  -  0)]  +  cos  [4 irf0t  -  2irf0(a  +  0)  +  2 <j>]  jdt 

=  exp  [— i2x/(or  -  0)]  |cos  [2x/0(a  -  0)]  / 

f  sin  [4x/p<  —  2x/0(r  +  /?)  +  2<t>)  '(j+1)' 


4ff/o 


j/ 


}  U-8) 


which  if  /  l/4x/0 


!  exp  [-«2x/(a  -  /?)]  cos  [2x/o(cr  -  0)]  I 


(1-9) 


Now  we  should  note  that 


exp  [—  i2vf(a  —  0)]  cos  [2x/o(a  -  0)\  I 

=  ^exp  [~*2x/(a  -  j3)]  {exp  [»27r/0(a  -  0)]  +  exp  (-i2irJQ(a  -  /?)]} 

=  j{exP  [-*2x(a  -  0)(f  -  f0)]  +  exp  [-t27r(a  -  0)(f  +  /o)]}  (I  -  10) 


Substituting  (1-10)  in  (1-7)  we  get 


Xjt  «  —r~  j  [  [  H(f  -  fk)w(a)w*(0){exp[-i2x(a  -  0)(f  -  f0)\ 

’  «/  — OO  J  —  oo  J —oo 

+  exp[-i2x(a  -  0)(f  +  f0)]}da  d0  df  (I  -  11) 


If  we  define 


G(f)  = 


r 

I  w(z)exp(—i2itfz)dz 

J  —  OO 


(I  —  12) 


which  is  the  magnitude  squared  of  the  Fourier  transformed  window  function,  then  we  can  rewrite 


G(f)  =  f  w(a)exp(-i2irfa)da  f  w‘(0)exp(i2irf0)d0 

J —OO  J— OO 

/OO  fOO 

•oo  j  —  oo 


w(a)u>‘(0)  exp(-i2nfa )  exp(i2ir  f  0)da  d0 
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and  hence 


and  similarly 


G(f  +  fo)=  f  f  w(a)w*(0)exp[-i2ir(a- /?)(/  +  f0)]dad{3  (1-14) 

J  —  oo  J  —  oo 


/oo  r  oo 

/  w(ot)w*  (P)  exp  [— t'27r(a  —  /?)(/  —  /0)]  da  d/3  (I-  15) 

’OO  J  —  OO 


Substituting  (1-14)  and  (1-15)  in  (1-11)  we  get 


A2 1  r°° 

Xik  *  —■  /  //(/  -  /t)[G(/  -  /o)  +  G(f  +  /0)]d/ 

^  J  —  OO 


Now  if  we  assume  that  //(/)  is  symmetrical  about  /  =  0,  then 


/OO  rCO 

W(/  -  fk)G(f  +  /0)d/  =  /  H(fk  -  f)G(f  +  f0)df 

-oo  •/  — OO 

r°° 

=  /  tf[/*-(/  +  /o)  +  /o]C(/  +  /o)d/ 

*/ -OO 

=  r  H(fk  +  /o  -  f')G(f')df'  =  H(fk  +  /o) 

J  —  OO 

where, 

r°° 

n(f)=  /  H(f  —  f')G(f  )df 

J  —  OO 

whicli  is  the  convolution  between  the  functions  G  and  77.  And  similarly,  it  can  be  shown  that 

«(/*  -  /o)  =  /°°  //(/  -  fk)G(f  -  fo)df 

J  —  OO 

Hence,  substituting  (1-17)  and  (1-19)  in  (1-16)  we  get  that 


Xjk  «  — —  {‘H-Uk  -  fo)  +  "H{fk  +  /o)} 

But  ‘H(fk  4-  fo)  is  very  small  compared  to  ~H(fk  —  fo)  for  large  values  of  fk  or  fo,  hence 


A2  / 

“  fo) 


(1  -  16) 


(7-17) 


(7-18) 


(7-19) 


(7  —  20) 


(7  —  21) 
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From  equation  (5)  we  have  that 


w(z )  =  exp  —az  —  4 T2  ^  rect  ^ 


(II-  1) 


or  we  could  normalize  w(z )  and  get 


w(x)  =  exp  —Lx  —  4 T2  ^x  —  ^  rect  ^ 


(II  —  2) 


where  L  =  ar  and  x  =  £ . 


From  equation  (4)  and  (II-l)  we  have  that 


G(f)  =  |y  tn(z)exp(— i2nfz)dz 


(II  -  3) 


A  closed  form  solution  to  this  equation  does  not  exist  in  general  unless  some  simplification  or  approx¬ 
imation  is  done.  The  rectangular  rule  can  provide  a  numerical  estimate  of  this  integration.  This  gives  the 
following  equation  as  an  estimate  of  equation  (II-3): 


Ga(f)=  A t  w(nAt)  exp(— »'2;r/nA<) 


(II  —  4) 


Now,  in  this  case  we  will  let  At  =  r/N  where  N  is  the  number  of  points  from  tv(z)  that  will  be  used  to 
calculate  one  sample  of  G(f). 

Replacing  At  =  r/N  in  (II-4)  we  have: 


n  — n  '  ' 


(II-  5) 


if  we  let  v  =  fr  we  get 


„  ,  ,  r  1  /nr\  /  —  t2nvn\ 

c.w=  j?  E  “>  (-) (— —) 

n  =  0  x  ' 


(II -6) 


or  if  we  express  the  above  equation  in  rectangular  coordinates,  we  get 


Ga(v)  =  \tI<\2 


(II-  7) 


where, 


r  EnJo  h  M  cos  (JT1  -  iw  ill sin  (fgOl 

N 


(II-  8) 


hence, 


G.(v)  =  r2  {[^(/v)]2  +  [/m(/C)]2} 


To  consider  the  error  associated  with  the  approximation  of  (II-9)  let 


N-l 


Xa(f)  =  At  ^2  u>(nAt)exp(—i2wfnAt) 


(7/ -10) 


n=0 


It  is  easy  to  see  from  (II-4)  that 


Ga(f)  =  |Xa(/)|2 


(77-11) 


We  note  that  the  Fourier  transform  of  a  delayed  Dirac  delta  function,  S(t  -  <0),  is  e~,2*Jto.  Thus  we 


may  regard  (11-10)  as  the  Fourier  transform  of  At£„J0l  w(nAt)6(t  -  nAt): 


Xa(f)  =  T  | At  w(nAt)6(t  -  nA<)  j 


(77-12) 


Since  tv(z)  is  zero  outside  the  interval  [0,  r],  we  may  rewrite  (11-12)  as 


Xa(f)  =  T  |  w(z)  ■  At  jr  «(z-nA<)l 

l  n=  — oo  ) 


(77-13) 


which  can  be  transformed,  using  the  well  known  property  for  the  Fourier  transform  of  multiplied  time 
functions,  to  the  following: 


»(/)  =  w(/)^|a<  y.  *(*-»**)} 

l  nsr-oo  ) 


which  becomes 


*.</>  =  ww  £«(/-£) 


which  finally  becomes 


*.</)= 

ns-oo 


(II  —  14) 


So  we  see  that  the  approximation  of  equation  (11-9)  is  in  fact  an  aliased  version  of  G(f).  We  also  see 
that  Ga(f)  is  periodic  with  period  1/A t,  hence  there  is  no  point  in  calculating  G„(f)  outside  the  range  of 
frequencies 


-N  . 
—  <f< 


N_ 

It 


since  A t  =  r/N .  Or  if  we  use  the  normalized  version  Ga(v),  then  the  range  of  frequencies  is 


(tv  7«.  c.«tificr 
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